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ABSTRACT 

Neutron-star cores may be hosts of a unique mixture of a neutron superfluid and a proton superconductor. 
Compelling theoretical arguments have been presented over the years that if the proton superconductor is of 
type II, than the superconductor fluxtubes and superfluid vortices should be strongly coupled and hence the 
vortices should be pinned to the proton-electron plasma in the core. We explore the effect of this pinning 
on the hydromagnetic waves in the core, and discuss 2 astrophysical applications of our results: 1. We 
show that even in the case of strong pinning, the core Alfven waves thought to be responsible for the low- 
frequency magnetar quasi-periodic oscillations (QPO) are not significantly mass-loaded by the neutrons. The 
decoupling of ~ 0.95 of the core mass from the Alfven waves is in fact required in order to explain the QPO 
frequencies, for simple magnetic geometries and for magnetic fields not greater than 10 15 Gauss. 2. We show 
that in the case of strong vortex pinning, hydromagnetic stresses exert stabilizing influence on the Glaberson 
instability, which has recently been proposed as a potential source of superfluid turbulence in neutron stars. 



Subject headings: Neutron stars 



Introduction 



Since the late 1950's, it has been realized that neutron- 
star interior may consist of a number of quantum fluids 
(see Shapiro & Teukolsky 1983 for a review). Currently, 
it is thought that both neutron superfluid and proton su- 
perconductor are likely to coexist in the neutron-star cores 
(see, e.g., Link 2007 for a discussion). Several researchers 
have argued that if the proton superconductivity were of 
the type II, then the superconductor fluxtubes would cou- 
ple strongly to the neutron superfluid vortices. This line of 
reasoning is based on the fact that nuclear forces contain 
velocity-dependent terms, which results in the entrainment 
of protons in the neutron supercurrent (Alpar, Langer, 
& Sauls, 1984). Therefore, the vortices are sheathed by 
charged currents entrained in the superfluid flow, and are 
strongly magnetized. Magnetic fluxtubes interact strongly 
with the magnetized vortices, similar to the way in which 
the fluxtubes interact between each other (Ruderman, Zhu, 
& Chen 1998, and references therein). As a result of this 
coupling, the vortices get strongly pinned to the proton- 
electron plasma in the core. Such pinning would have 
important implications for the neutron-star phenomenol- 



ogy. Ruderman, Zhu, & Chen (1998) have argued that the 
vortex-pinning in the core may be responsible for the ob- 
served glitches in the pulsar rotation rates. Link (2003) 
has considered the effect of the vortex-fluxtube interac- 
tion on the dynamics of the precessing pulsar PSR 1828- 
11 (observed by Stairs, Lyne, & Shemar 2000). Building 
on the theoretical work by Shaham (1977) and Sedrakian, 
Wasserman, & Cordes (1999), he has concluded that the 
interaction, if present, would ultimately lead to the fast 
precession. Since PSR 1828-11 is precessing slowly and 
persistently, Link (2003) has argued that the core vortex 
pinning is excluded by the observations and hence that ei- 
ther the proton superconductor might be of type I, or that 
both proton and neutron condensates do not coexist inside 
that pulsar. While Link's argument is suggestive, we be- 
lieve it is premature to rule out strong vortex pinning in 
the cores of all neutron stars. 

In this paper we consider hydromagnetic waves in the 
case when the neutron vortices are strongly pinned to the 
proton-electron plasma in the core. We have 2 main astro- 
physical motivations for studying this problem. 

The first one is due to the fairly recent observations 
of the quasi-periodic oscillations (QPOs) of the x-ray lu- 
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minosity in the tails of giant magnetar flares (Israel et 
al. 2005, Strohmayer & Watts 2005, 2006, see also ear- 
lier but lower signal-to-noise measurements of Barat et 
al. 1983). Israel et al. 2005 has argued that the lowest- 
frequency and the longest-lived QPO of ~ 18Hz is likely to 
represent an Alfven-type oscillation in the magnetar core 
(this frequency is too low to be associated with the tor- 
sional modes of the crust). Levin (2006) has shown that 
for a magnetar-strength field the crustal motion [which is 
thought to be either powering power the flare (Thompson 
& Duncan 1995) or responding to a global reconnection 
event in the magnetosphere (Lyutikov 2003)] would excite 
the core Alfven waves on the timescale of several oscilla- 
tion periods. Since then, a significant body of theoretical 
work has been devoted to a study of global magnetar vi- 
brations, which would involve both hydromagnetic waves 
in the core and elasto-magnetic shear waves in the crust 
[Glampedakis, Samuelsson, & Andersson 2006, Levin 2007 
(from here on L07), Sotani, Kokkotas, & Stergioulas 2007, 
Lee 2007]. In particular, L07 has argued that the long- 
lived low-frequency QPOs are associated with the special 
spectral points of the Afven continuum in the magnetar 
core. For simple B-field configurations these special points 
can be found analytically, and do not depend on the details 
of the crust. For example, for a uniform internal B-field, 
the lowest QPO is expected at the frequency 
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where R is the radius of the fluid part of the star, p c is the 
density of the the part of the fluid which is coupled to the 
Alfven waves, and B c g is the effective magnetic field which 
is given bjQ 

B c s = V BB CT . (2) 

Here B is the average magnetic field, B cl - ~ 10 15 G is the 
value of the critical magnetic field confined to the flux- 
tubes. ^From Eq. (f5]) we see that the magnetar QPO 
frequencies could provide an interesting constraint on the 
magnetic-field strength and geometry. However, interac- 
tion between neutron and proton superfluids could affect 
core Alfven waves, by effectively mass-loading them with 
neutrons. We shall consider the extreme case of such 



lr The occurence of B c g in Eq. and |21l can be understood as follows: 
the magnetic tension force acting on surface S perpendicular to the 
fluxtubes containing N fluxtubes is given by NAS- B^ t /4tt = £■ T e g , 
where AS is the cross section of a single fluxtube, and T e fj is the 
effective tensile stress. The magnetic flux through E is $ = S- B = 
ATAS-Bcr, and we find T cff = BB CI /4n = B^ s /4n (as opposed to 
T = B 2 /47r, where T is the corresponding part of the Maxwell stress 
tensor). (A detailed derivation is given in Easson & Pethick, (1977)). 



interaction — the strong vortex pinning on the fluxtubes, 
and show that it does not significantly alter the Alfven- 
wave propagation in slowly-spinning magnetars (but is im- 
portant for the Alfven waves in the fast-spinning radio- 
pulsars). This simplifies the interpretation of the QPO 
frequencies, and shows that it is valid to assume that p c 
is just the density of protons, about 5% of the total fluid 
density. 

Our second motivation is the recent theoretical discus- 
sion of the superfluid turbulence in the neutron-star cores. 
Superfluid turbulence has been discussed in the context of 
the laboratory Helium fluid for the last 30 years (see, e.g., 
Donnelly 1991 for a review). In a ground-breaking series 
Peralta, Melatos, Giacobello, & Ooi at the University of 
Melbourne (2005, 2006; hereafter PMG05 and PMG06) 
have applied the superfluid-helium ideas to neutron-star 
interiors. PMGO have developed from scratch a code 
which studies numerically the 2-component superfluid dy- 
namics. The 2 components in PMGO are the neutron su- 
perfluid and the normal neutron fluid which are coupled 
via the mutual friction force at the superfluid vortices; this 
mixture is expected if the core temperature is a signifi- 
cant fraction of the critical temperature of the superfluid. 
The equations of motions used by PMGO were derived 
by Hall and Vinen (1956) and Bekharevich and Khalat- 
nikov (1961). PMG05 have studied, for the first time, the 
superfluid spherical Taylor-Couette flow, and find that it 
becomes turbulent in 2 steps: 1. the normal component 
develops meridional circulation due to the Eckman pump- 
ing (see, e.g., Pedlosky 1987), and 2. the component of 
the meridional flow of the normal fluid which is directed 
along the superfluid vortices drives the vortex Kelvin waves 
unstable; this is known as the Glaberson (or sometimes 
Donnclly-Glaberson) instability (Glaberson, Johnson, Os- 
termeier 1974, Donnelly 1991). In PMG05 simulations, 
the Glaberson instability leads to turbulence. Interest- 
ingly, PMG06 and Melatos and Peralta (2007) demon- 
strate that the superfluid turbulence could affect the pul- 
sar spin and could be behind the well-known pulsar timing 
noise and spin glitches. 

More recently, Sidery, Andersson, and Gomer (2007, 
hereafter SAG), and Glampedakis, Andersson, and Jones 
(2007, hereafter GAJ1, and 2008, hereafter GAJ2) devel- 
oped an analytical theory of the Glaberson instability in 
neutron stars. Their 2-component fluid consists of the neu- 
tron superfluid and the proton superconductor, which are, 
like in PMGO, coupled via the mutual- friction force. SAG 
have considered the limit of the weak mutual friction (see 
below) and infinitely heavy proton superfluid, and found 
that it was the inertial waves in the neutron superfluid 
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which were subject to the Glaberson instability. GAJ1 and 
GAJ2 have extended this analysis to include the regime 
of realistic proton-to-neutron mass ratio and of arbitrar- 
ily strong mutual friction. Notably, they found that the 
Glaberson instability operated even in the regime of strong 
pinning. But where would the initial relative flow between 
the protons and neutrons come from? GAJ have argued 
that if the pinned neutron vortices were misaligned with 
the angular velocity of the protons, then this would natu- 
rally lead to the relative proton-neutron flow which would 
potentially be strong enough to drive the turbulence in 
some parts of the star. Without turbulence, the star with 
misaligned pinned vortices would undergo fast precession 
(Shaham 1977) which, although probably hard to detect, 
has yet not been observed in any of the currently-timed 
radio pulsars or magnetars. GAJ2 argue that this the 
Glaberson-instability-induced turbulence may generically 
prevent the star getting into a configuration with the fast 
precession. However, this conclusion is premature. One 
important piece of physics which is not considered in GAJ 
is the strong hydromagnetic stress inside the proton super- 
fluid, which, as we show below, has a suppressing effect on 
the Glaberson instability and hence on the development 
of the superfluid turbulence. We will derive, however, a 
robust upper limit on the angle of fast precession, which 
is determined by the maximum possible mutual torque be- 
tween the neutron superfluid and the proton superconduc- 
tor. The maximal precession angle turns out to be much 
smaller than 1 degree, and thus it is not surprising that 
the fast precession has never been observed in real neu- 
tron stars. 

The plan of the paper is as follows. In the next section 
we derive the dispersion relation for the hydromagnetic 
waves when the superfluid vortices are rigidly attached to 
the core plasma. In sections 3 and 4, we consider applica- 
tions to oscillating magnetars and precessing pulsars, re- 
spectively. We end with the general discussion in section 
5. 

2. Dispersion relations 

As a starting point, we utilize the plane- wave analysis 
outlined in GAJ. We follow closely the notation of and 
reasoning behind GAJl's basic equations (1) — (7), and our 
derived dispersion relation is identical to their Eq. (10) in 
the limit of zero hydromagnetic stress, but has important 
extra terms when the hydromagnetic stress is included. We 
begin with the two-fluid dynamical equations, cf. Eqs. (1) 



and (2) in GAJ1: 

D?v n + V^„ = 2v n x (1 + f nii (3) 
D^v p + Vip p = 2v p x - f m f/xp + v ce V 2 Vp 

+ /hm (4) 

Here v n and v p are the velocities of the neutron and proton 
condensates, respectively, the full time derivatives D t are 
defined in the usual way as — d/dt + v, ltP ■ V, tpn.p is 
the sum of specific chemical and gravitational potentials, 
/ m f is the acceleration of the neutron superfluid due to 
the mutual friction between its vortices and the charged 
plasma, x p = p v j p n is the density ratio between the proton 
charged and neutral components of the interior (~ 5%), f2 
is the angular velocity of the rotating frame in which all of 
the velocities are defined, v ee is the kinetic viscosity of the 
plasma due to electron-electron scattering, and 

Am = B cS - VB cS /4irpp (5) 

is the acceleration of the plasma due to the hydromagnetic 
stress. We note that because in a type-II superconductor 
the distance between the fluxtubes is much larger than the 
fluxtube diameter, we ignore magnetic pressure. In writ- 
ing down Eq. ©, we have followed GAJ and neglected 
explicitly the effect of superfluid entrainmcnt (which we 
expect will not qualitatively change our results) and the 
individual tension force for the vortices (which is can be 
neglected if the wavelength of the waves is much greater 
than the inter- vortex distance). We shall use the follow- 
ing conventional form (Hall & Vinen 1956, Bekharevich & 
Khalatnikov 1961, PMGO, SAG, and GAJ) for the mutual- 
friction force: 

-> R R 2 

/mf = j — R2 &n X (u„ X UJ np ) + 1 + R2 ^ n X "kpj ( 6 ) 

where (i„ = Vx V n is the vorticity of the neutron fluid in a 
non-rotating frame (here V n is the neutron velocity in the 
non-rotating frame jj, ui n — w n /|cJ^| is the associated unit 
vector, w np = v n — v p , and R is the dimensionless number 
measuring the strength of the drag between the neutron 
vortices and the plasma. When R -C 1 (the weak-drag 
limit), the first term on the right-hand side dominates. 
This entails that the neutron vortices mostly follow the 
motion of the neutron superfluid in the direction perpen- 
dicular to uo n . When R ^> 1 (the strong-drag limit), the 
second term on the right-hand side dominates. This en- 
tails that the neutron vortices mostly follow the plasma 



2 In GAJl, uj n is erroneously defined as Vx« B . However, in their 
subsequent calculations they, most likely, use the correct expression. 
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motion. When R — oo, which is the case on which this 
paper focuses, the vortices get pinned to the plasma. In 
this limit, the plasma and the neutron superfluid interact 
exclusively via the Magnus force arising from the relative 
motion between the neutron vortices and neutron super- 
fluid. 

We choose the background state as follows: 1. the z-axis 
is directed along 17; 2. the neutron vortices are aligned with 
51 = 51e z , and are at rest in the rotating frame; 3. in the 
same frame, the plasma has a background velocity wq = 
woe z , which is directed along the vortices; 4. the mean 
magnetic field is directed along the vortices, B — Be z . We 
consider waves which are propagating along the z-axis. We 
are interested in the waves for which the restoring force 
is the combination of hydromagnetic stress, the Coriolis 
force, and the Magnus force. This means that the wave 
must be nearly incompressible, which implies 



k ■ Sv n p = 0. 



(7) 



Here k is the wavevector, 8v nyP is the neutron/proton ve- 
locity perturbation due to the wave. Incompressibility 
and assumed homogeneity of the background state imply 

Slpn.p = 0. 

Let us introduce the Lagrangian displacement vectors 
£n,p of the neutron and proton fluids from their background 
positions, with 8v n>p = -D™' p £ njP . We are looking for the 
solutions of the form 

£n, P {z,t) = [Co P e x +Q P e v ] e^ t+kz \ (8) 

where a is the angular frequency of the wave. 

We now perturb Equations 10]), and (0; we set 
R = oo in the latter. To the linear order in the velocity 
perturbation, we have: 

DfSv n = d 2 C/dt 2 = -a 2 £„, 
D?8v p = -(<T + kw ) 2 £ p , (9) 
6f m f = 2Q x ^6v n - Svj^J - (V x 8v n ^ x w Q , 

v ee S7 2 8v p = -iv cc k 2 (<7 + kw )£ p , 
6 flm = c 2 A d%jdz 2 = -c\k%. 



Here ca = \J BB CI / (47rp p ) is the Alfven velocity in the 
plasma. The expression for 5fh m is obtained using the 
magnetic induction equation. Substituting these into 
Eqs. and (0]), and using V x £ = ik x £ together with 
Sv n = ia^ n and Sv p — i(a + fcu>o)£ p , we get 2 linear vector 
equations for and £ p . It is now convenient to proceed 
as follows: 



Let us represent a vector £ = ^e x + (,y£y by a complex 
number £ = ^ x + . Then a vector e z x £ is represented 
by i£. By using this, we can immediately rewrite the 2 real 
vector equations as 2 complex scalar equations: 



rr + 2<> ( i - — ) - (iVee + c\/a) k 



257 — kwo 



0, 

(10) 



cr£„ = 0, 



where a — a + kwo . This pair of equation yields immedi- 
ately the complex dispersion relation: 



251 1 



i 

Xp 

251(251 - kw ) 

x p 



2 1.2 



c A k 



0. (11) 



The dispersion relation for arbitrary R is derived, for com- 
pleteness, in Appendix A. 

In the next 2 sections we consider 2 applications of the 
relation Eq. (fTTjl . 

3. Hydromagnetic waves in magnetars 

In this section we assume that there is no 51-directed 
relative proton-neutron flow, i.e. we assume wq = 0. We 
also set v ee to zero, since the ratio of the viscous to hydro- 
magnetic stress is given by 



v cc a/c 2 A < 1. 



(12) 



With these simplifications, the dispersion relation (fTT|) 
gives 



= -£11 




It is important to note that in this expression ca is a func- 
tion of only the proton density p p (c\ = BB cr /A-Kp p ). All 
observed magnetars are slowly rotating, with 51 ~ lrad/s. 
The observed lowest angular frequency for a magnetar 
QPO is 18Hz, thus a ~ 113rad/s. The sum of Magnus 
and Coriolis forces, represented by the terms with 51, con- 
tribute only a fraction 8a /a to the wave frequency, given 
by 



8*/*,^- = 0.1 8 ( " ) 

Xptr \ lrad/s J \ a 



0.05 x 



(14) 
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We note that this does not depend on the assumption that 
a represents some fundamental Alfven mode. 

^From the above equation, it is clear that for hydro- 
magnetic waves associated with the magnetar QPO fre- 
quencies of 18Hz and higher, the magnus forces from neu- 
tron superfiuid introduce only a small correction to their 
propagation. Thus we conclude that the magnetar oscilla- 
tions (as seen in the giant-flare QPOs) even in the case of 
strong pinning, do not couple significantly to the neutron 
superfluidd Therefore, given the knowledge of the inter- 
nal magnetic field, one should use p c ~ x p p n in Eq. (fT|) to 
determine the frequency of the lowest QPO which, accord- 
ing to Levin (2007), corresponds to the turning point of 
the Alfven continuum in the core. For Levin's the simplest 
computable magnetar model (uniform internal magnetic 
field and density), with the typical magnetar parameters, 
B = 10 15 G, R = 10km, p = 10 15 g/cm 3 and x p = 0.05, 
Eq. (fT]) gives v a ~ 22Hz, which is in qualitative agreement 
with the observed 18Hz. If the whole neutron superfiuid 
would mass-load the wave, this frequency would be reduced 
by a factor of ~ 4. While suggestive, the numbers above 
certainly should not be taken as evidence of neutron su- 
perfluidity, since the strength and geometry of magnetic 
fields inside the magnetar are highly uncertain. 

4. Application to precessing neutron stars 

Consider now a precessing neutron star where the neu- 
tron angular velocity and the crust+proton angular veloc- 
itjH are equal in magnitude O but are misaligned by an 
angle 9. Suppose that this angle is fixed due to the strong 
vortex pinning. The relative velocity of the proton super- 
fluid along the vortices is given by 



wq = £1 sin 9x2, 



(15) 



where X2 is the coordinate measured along VL n x fl p . Note 
that this expression agrees with Eq. (18) in GAJ1 when 
one notes that for small 9 their wobble angle 9 W equals 
I p 9 / I ni where I p and /„ are the moments of inertia of the 
proton and neutron components, respectively. 



3 We note that we have used the plane- wave analysis for what is likely a 
non-plane- wave oscillation. This is clearly a limitation of our formal- 
ism. However, the plane-wave analysis illustrates the physics which 
is also valid for oscillatory mode of any geometry, namely that for 
high-frequency waves Alfven restoring forces are much greater than 
the Magnus forces. This occurs essentially because the Magnus force 
is proportional to the velocity and thus scales as cr, while the to- 
tal restoring force scales as cr 2 . Thus, we believe that our analysis 
is qualitatively correct in the high-frequency regime, for non-plane- 
wave Alfven-type oscillations. 

4 The crust and the core protons are co-precessing; this is enforced by 
the hydro-magnetic stresses (Levin & D'Angelo, 2004). 



4.1. Glaberson instability criterion 

It is convenient to express the general solution of 
Eq. ((TTJ) as follows: 



a = —ku)Q —Oil 



iv v . v k 2 



± 



(16) 



D, 



where 

D = O 2 I 1 



+ c\k 2 



2kw n v 2 k A 



~k 2 n 1 1 - — 



1 



(17) 



This is essentially the same as Eq. (10) of GAJ1 when 
c A = 0. 

First, lets consider the non-viscous case with v cc = 0. 
Then the minimum of D is attained for k = Qwo/(x p c A ) 
and equals 



Anin =n 2 1 



1 



xp\c A 



Thus the Glaberson instability appears only for 
w > ca{1 + x p ) ~ c A . 



(18) 



(19) 



We now allow for the viscous term in Eqs. (jTHJ) and (jTTJ) . 
We find that a weak, viscosity-driven instability appears 
at a smaller velocity 



(20) 




W ± \Jwq~ 4c A Xp 



see Appendix B for the mathematical details. 



(21) 
(22) 



Equation (|20|) expresses the lowest relative proton- 
neutron velocity which is required for the initiation of the 
Glaberson instability. Whether this velocity is acheived 
depends on the misalignment angle 9 between the proton 
and neutron anglular velocity vectors. In the next sub- 
section, we derive a simple but rigorous upper bound on 
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4.2. The maximum misalignment angle and rela- 
tive flow velocity for fast precession. 

The misalignment angle 9 can be constrained, by requir- 
ing that the precessional torque t p of the proton compo- 
nent not exceed the maximum torque r m that the vortices 
can exert on the fluxtubes through magnetic stresses. The 
precessional torque is given bjd 

T p = L n tt p sin# (23) 

Here L n = IQ is the proton angular momentum, / is the 
total stellar moment of inertia and tt p = f2. We find that 
for a typical neutron star with the mass of M — 1.4M Q 
and radius of R = 10 6 cm, the precessional torque is given 

by 

t p ~ 4 • 10 46 sin6>(P/ls)- 2 gcm 2 s- 2 . (24) 

Here P is the neutron-star rotational period. 

The maximal physically-admissible torque r m can be 
computed by assuming that the vortices have maxi- 
mal contact with the fluxtubes, i.e. that the vortex is 
pushed/pulled with the stress of B 2 r /(8n) accross its whole 
side surface. The maximal torque exerted on a single vor- 
tex is given by 

t v = |^Z 2 d, (25) 

where d ~ 2- 10~ 12 cm is the vortex diameter, and I is the 
vortex half-length. The total number of vortices is given 
by 

N = ttR 2 ti v ~ 3- 10 16 (.H/lQ 6 cm) 2 (P/ls) _ \ (26) 

where n v is the per-area vortex density (Shapiro and 
Teukolsky 1983, Link 2003). For a spherical star with a 
dense grid of the linear vortices, the average value of I 2 is 
P 2 /2. Thus, by combining Eqs. (25) and (26), we arrive 
to the following expression: 

Tm = -rfR 4 dn v ~ 1.3 • lO^P/ls^gcmV 2 . (27) 

^From Eqs. (|27|) and (|24|) . we see that our requirement 
T m > Tp implies that 

6 < 2°(P/ls), (28) 



Since the neutrons are pinned to the protons, the torque acting on the 
neutrons is given by the cross product of the instantaneous angular 
velocity of the protons and the neutrons angular momentum, and 
is therefore independent of x p . In our derivation we assume that 
the angular velocities of the protons and the neutrons have the same 
magnitude fi. 



and therefore 

w ~ flOR < 2- 10 5 cm/s. (29) 

This upper limit on wq is spin- independent. 

So is this velocity sufficient to drive the Glaberson in- 
stability? Equation (f20"|) tells us that for x p — 0.05, 
B = 10 12 G and p n = 10 15 g/cm 3 , the critical relative 
velocity is wq ~ 6 • 10 5 cm/s. Thus we conclude that in 
the presence of strong vortex pinning and magnetic fields 
B > 10 11 G the misalignment between the proton- and neu- 
tron angular velocities is unlikely to become large enough 
to provide wind velocities wq that are sufficient to drive 
the Glaberson instability. 

5. Outlook 

The calculations of this paper have 2 main astrophysical 
implications. First, we have shown that the Alfven waves 
which are associated with magnetar QPOs are not signifi- 
cantly mass-loaded by a neutron superfluid, even if the su- 
perfluid vortices are strongly pinned to the proton-electron 
plasma. For B — 10 15 G and the simplest B-field geome- 
try, the expected frequency of lowest magnetar QPO is in 
remarkable agreement with observations, if only the pro- 
tons, i.e. about 0.05 of the core mass, are mass-loading the 
Alfven waves. Strong vortex pinning will, however, have a 
strong effect on the Alfven waves in rapidly spinning and 
relatively non-magnetic neutron stars, i.e. those ones in the 
Low-Mass X-ray Binaries. In these stars the Alfven waves 
may play an important role in damping of the r-mode in- 
stability, as discussed by Mendell (2001) and Kinney and 
Mendell (2003) for the cases of non-superfluid and super- 
fluid core, respectively. Kinney and Mendell (2003) had 
not considered the vortex pinning (see also Mendell 1991); 
however from Eq. (fT4")) and from the fact that a ~ fi for an 
R-modc, it is clear that the strong vortex pinning would 
have a huge (of order l/x p ) effect on the Alfven waves with 
the R-mode frequency. 

Second, we have shown that the hydromagnetic stresses 
generally suppress the Glaberson instability in the proton- 
neutron superfluid mixture, in the case of strongly pinned 
vorticefl We have also provided a robust upper bound 
Eq. |28|) on the angle between proton and neutron angu- 
lar velocities in the fast-precessing neutron stars. Even for 
slowly-spinning magnetars, the misalignment angle cannot 
exceed 20 degrees, which implies a wobble angle no greater 
than 1 degree. Thus a detection of neutron-star fast pre- 
cession is difficult, if not impossible, observationally. An 



We have not considered here the PMGO case when some normal 
neutron component is present and is driving the instability. 
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inspection of the XMM timing data on known anomalous 
x-ray pulsars produces no statistically-significant periodic 
signal which could be interpreted as fast precession [van 
Kerkwijk 2008, private communications]. 
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A. Dispersion relation for arbitrary drag 

In this Appendix we perform a plane- wave analysis of the two- fluid dynamical equations defined in eq. (3) and (4), 
and derive a dispersion relation for arbitrary drag strength R. Using the plane wave solutions from equation (7), we can 
perturb equations (3) and (4). Retaining linear terms we get: 



R 



mf 



1 + R 2 



D?5v n = -a 2 l L 

DfSv p = -v 2 i p 
25v n x = 2iOer£ n x e z 
2Sv p X Q = 2i^la^ p x e z 

2iOe z x (e z x (a£ n - a^p^ + kw ae z x ((e z x £j x e^j 
R 2 



+ 



1 + R 2 



2ifle z x (a!; n - cr^ - w kae z x (e z x £ n ^J 



u ee \7 2 Sv p = -v cc k 2 a^p 



Sfhr 



(Al) 

(A2) 
(A3) 
(A4) 

(A5) 



(A6) 
(A7) 



Here a = a + wgk and ca — y BB cr / AiTp c is the Alfven velocity in the plasma. We can simplify these expressions by 
using the same trick as in Section 2: Let us represent the vector £ = £ x e x + £, y e y by a complex number £ = £ x + i^ y - 
The cross product e z x £ corresponding to a simple rotation in the xy-plane, can then be represented by i£. Using this, 
we convert the two real vector equations (3) and (4) into two complex scalar equations: 



2ClCa£ p = [a + (C-l) {w a k - 20)] tr|„ 



20 1 



C 



-A ,2 



c 



(wok — 20) a^ n 



(A8) 
(A9) 



where C = ^ ■ This pair of equations yields the following complex dispersion relation: 



^•3 , -2 

a + a 



( 20 
(w k - 20) (C - 1) (20 - iv cc k 2 ) + — 



20C\ 
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2 t,2 



c 2 A k 2 (C - 1) (w /s - 20) = 



In the strong drag limit i? — > oo (C =1) this cubic equation simplifies significantly 

"2O(2O-u> fc) 



a 2 - 




2 " (^) 

















= o, 



which is Eq. (fTTjl in the text. 



(A10) 



(All) 
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B. Instability criterion for non-zero viscosity 

In this Appendix we derive a criterion for instability in the case of non-negligible viscosity, i.e. Eq. (|20j) . We rewrite 
Eq. (fTB]) as follows: 



where 



a = A + iB± VC + iD (Bl) 

A = -w k - fl ( 1 - — 
V x p 

L> 



2 



C = n 2 (l + ±-)\c 2 A k 2 - 2 -^-B 2 (B2) 



D = -20 ( 1 1 13 

The state of marginal stability is given by 



Im (a) — B ± Im (ycTiD^j = 0. (B3) 



In order to find a convenient expression for Im (^JC + iD^j we write (yj C + iD^j in polar form. The imaginary part is 
then given by 

C 



Im[VC + W) = (C 2 + D 2 ysin ( - arccos ( — = ^ ) ) (B4) 



B 2 = ±VC 2 + D*-±C (B5) 



Combining this with Eq. IB3P , we find 



Taking the square of this expression, and using Eqs (|B2|) we arrive at 

c\k 2 = (B6) 



2<> (20 - /(•„/,•) i 2 , , 



And therefore, 

fe± = —^2" 

XpCj^ 



w ± \Jwl~ 4c 2 A Xp 



(B7) 



The unstable waves have wave- vectors in the interval fc_ < k < fc + , provided that fc_ and k + are real. Thus the 
criterion for the instability is 

w Q > 1ca^[x~p, (B8) 

this is Eq. (|20p of the paper. We now make an estimate of the growth rate in this instability window. For a realistic 
neutron star, we take v ee ~ 10 9 cm 2 s~ 1 at T ~ 10 7 K (Flowers & Itoh 1979, Cutler & Lindblom 1987, Andersson, Comer 
& Glampedakis, 2005), rj 27rrad/s, ca ~ 10 6 cm/s and therefore for wo ~ ca, the wave-vector of an unstable wave is 
k ~ 10 _4 cm _1 . We now note that u ee k 2 /2 « Q/x p . Therefore the terms B from Eq. (|B1[) and B 2 from Eq. (|B2|) have 
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a negligible contribution to Im (tr) and can be ignored. 



First the case where C < 0. We note that Im (a) is completely dominated by C so that we arrive at the criterion of 
Eq. (fT§|) again. With a growth rate of 



- Im (a) = VC (B9) 
Next consider C > 0, i.e. w$ < (1 + x p ) ca- By means of a simple analysis in the complex plane one can show that 



D T f \ D 

Y< -Im{a)< S J- 



(B10) 



For wo > there is a range of k where the instability occurs; see Eq. (f2TT) . Substituting k = k + into Eq. (|B10j) 

for the maximum k in this range, we find the growth rate of the instability 



Im (a) 



^ x p c A 



w + \/wg - 4c/ A x p 



(Bll) 



We note that because the viscosity is small, this growth rate is much smaller than the one in Eq. (|B9 
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